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Abstract 

In a finite temperature Thomas-Fermi theory with realistic nuclear interac- 
tions, we construct caloric curves for finite nuclei enclosed in a sphere of about 
4 — 8 times the normal nuclear volume. The specific heat capacity C v shows 
a peaked structure that is possibly indicative of a liquid-gas phase transition 
in finite nuclear systems. 
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The equation of state (E.O.S.) of nuclear matter with realistic effective interactions 
shows a typical Van der Waals type behavior and a critical temperature of ~ 15 — 20 MeV 
Supported by the experimental observation of a power law behavior in the mass 
or charge distribution in proton |4],|] and heavy ion induced reactions ||Q], ^he idea °f 
liquid-gas phase transition in nuclear matter or finite nuclear systems P,|5|-|T(J[] has gotten 
considerable interest in the literature. Theoretical speculations and possible experimental 
indications of a limiting temperature |TTHl5l in finite nuclei at ~ 5 — 7 MeV, above which 



the nucleus becomes unstable and breaks up into many fragments, also calls for a possible 
connection between the limiting temperature and the phase transition. Phase transitions 
are normally signalled by peaks in the specific heat at constant volume, C v as temperature 
increases. Fragmentation calculations in the microcanonical algorithm of Gross [16[] and in 



the Copenhagen canonical description [17,18] show such peaks. Recent calculations by Das 



Gupta et. al. |T9| in the lattice gas model for fragmentation also show such a structure. 
Renewed interest in this subject was further fueled by the recent experimental observation 
2T| in the caloric curve of a near constancy of temperature in the excitation energy range 



of ~ 4 — 10 MeV/nucleon in Au + Au collisions. This prompted us to find out whether the 
trends in the caloric curve as seen in the experiment or in fragmentation calculations are 
reproduced in a finite temperature Thomas-Fermi (TF) theory. To our knowledge this is the 
first calculation of its kind with a realistic effective interaction. In the context of an exactly 



solvable Fermion model, Rossignoli et al [21] have earlier calculated the specific heat of a 



finite nucleus in the grand canonical mean field theory with Lipkin's model hamiltonian, 
but found no structure in it as a function of temperature. The structure appeared in the 
canonical calculation, with inclusion of correlations. 

In our refined Thomas-Fermi (TF) model, the interaction density is calculated with a 
Seyler-Blanchard type |EJ momentum and density dependent finite range two-body effective 



interaction [13]. The interaction is given by 



Ves(r,p,p) = C hu [vi{r,p) +v 2 (r,p)} (1) 
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Wl = -(l-p 2 /& 2 )/(ri,r 2 ) 
V2 = d 2 [ Pl {r l )+p 2 {r 2 )] n f{v ll v 2 ) (2) 

with 

p-|n-r a | j n 

I ri - r 2 | I a 

Here a is the spatial range and b the strength of repulsion in the momentum dependence 
of the interaction, r =| ri — r 2 | and p = | pi — p 2 | are the relative distance and relative 
momenta of the two interacting nucleons. The subscripts I and u in the strength C refer to 
like pair (n-n or p-p) or unlike pair (n-p) interaction respectively, d and n are measures of 
the strength of the density dependence of the interaction and pi and p 2 are the densities at 
the sites of the two nucleons. 

The potential parameters are determined for a fixed value of n from a fit of the well- 
established bulk nuclear properties and the value of n is determined |T3| from a fit of the 
Giant Monopole Resonance energies over a broad mass spectrum. 

The Coulomb interaction energy density is given by the sum of the direct and exchange 
terms. They are given by 

e D (r) = e 2 Tcp p (r) J dr r' 2 p p (r')g(r,r'), (4) 

and 

e ex (r) = _^ (3 ^)^(r). (5) 
Here p p (r) is the proton density and 



/ f\ (r + r')- \ r-r'\ 

;(r,r) = . (6) 



With the potential chosen, the total energy density at a temperature T is then written as 

e{r) = p T (r)[T J 3/2 (rj T (r))/J 1/2 (rj T (r))(l - m* T {r)V^r)) + ~V T °(r)} (7) 

Here r refers to neutron or proton, the J's are the usual Fermi integrals, V® is the single 
particle potential ( for protons, it includes the Coulomb term), is the potential term 



that comes with momentum dependence and is associated with the effective mass m*. The 
fugacity r] T (r) is defined as 

rj r (r) = [ f i T -V T °(r)-V T 2 (r)]/T (8) 

where fi T is the chemical potential and V? is the rearrangement potential that appears for 
a density- dependent interaction. The total energy per particle at any temperature is then 
given by 

E(T) = J e(r) d 3 r/A. (9) 

Once the interaction energy density is known, the nuclear density can be obtained self- 
consistently and other observables of physical interest calculated. For details on the finite 
temperature TF theory, we refer to Ref. |13| . 



Since the continuum states of a nucleus at nonzero temperature are occupied with a 



finite probability given by a Fermi factor [p3|j , the particle density does not vanish at large 



distances. The observables then depend on the size of the box in which the calculations are 
performed. Guided by the practice that many calculations for heavy ion collisions are done 
by imposing that thermalisation occurs in a freeze-out volume, we fix a volume and find out 
the excitation energy as a function of temperature which allows for the determination of the 
specific heat at constant volume. 

We choose two systems, namely 150 Sm and 85 Kr. In the context of very heavy ion col- 
lisions at intermediate or higher energies, this mass range is of experimental interest. The 
calculations have been done for two confinement volumes, one at V = 4.0Vq and the other 
at V = 8.OV0, where Vq is the normal volume of the nucleus at zero temperature. The 
calculations at zero temperature are independent of the volumes taken; at low temperature 
of ~ 1 — 2 MeV, the observables are nearly independent of the volume. As the temperature 
increases, the central density is depleted. In Figure 1, the proton densities for 150 Sm cal- 
culated in the volume V = 8.OV0 are displayed for four temperatures, T = 5 MeV (dashed 
curve), T = 9 MeV (dotted curve), T = 9.5MeV (dash-dotted curve) and T = 10 MeV (full 
curve). At T = 5 MeV, the central density is depleted by ~ 4% compared to zero tempera- 



ture density, but has a long thin tail spread to the boundary. The behaviours at T = 9 and 
9.5 MeV are qualitatively the same, but with further depletion in the central density and 
a thicker tail. Beyond T = 9.5MeV, the change in the density starts being abrupt and the 
whole system looks like a uniform distribution of matter inside the volume. This is shown by 
a representative density distribution at T = lOMeV. The slight bump seen in the outer edge 
of the density is due to the Coulomb force. In Figure 2, the proton density at T = 10 MeV 
for the system at V — 8.0V (dashed curve) is compared with that calculated at V — 4.0Vo 
(full curve). The density calculated in smaller volume still shows a structure and the central 
density is depleted by only about 20% even at this high temperature. 

The excitation energy per particle E* is defined as E* = E{T) — E(T — 0). In Figure 
3, we display the caloric curve for the system 150 Sm. The upper dashed curve corresponds 
to V = 4.0Vo while the lower full curve corresponds to V = 8.OV0. At lower density, the 
excitation energy rises faster. For both volumes, initially the temperature rises faster with 
excitation energy, then its rise is slower. For the lower density, a kink is observed in the 
caloric curve at T ^ 10 MeV, after which the excitation energy rises almost linearly with 
temperature. For the higher density, the kink is much smaller and appears at a somewhat 
higher temperature. In Figure 4, the corresponding specific heats C v defined as 

C v = (dE*/dT) v (10) 

are displayed. Since we use units of MeV for both energy and temperature, the calculated 
C v is dimensionless. For both volumes, the specific heat shows a peak, the peak being much 
sharper for the case of a larger volume. For the smaller volume, the peak is at T ps 10.5 MeV 
while for the large volume the peak is shifted down by ps 1 MeV. We believe that the kink 
in the caloric curve or the peak in the specific heat are related to a phase transition in finite 
nuclei. From our calculations, we find that this transition temperature is weakly dependent 
on the confinement volume beyond V = 8V , e.g., for V as high as 20V , the transition 
temperature is shifted down further by only ps lMeV. The classical value of C v = 3/2 is 
reached at T ps 11 MeV for the case with V = 8.0V while for the smaller volume, it is 
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reached at T m 13 MeV. This is expected as the interaction becomes weaker either with 
increased volume or with increased temperature. 

In Figure 5, the caloric curve for the lower mass system 85 Kr is shown. The trends 
are nearly the same as in Figure 3. Figure 6 displays the specific heat for this system. 
In the calculation with V = 4Vo, a broad bump in the specific heat at T ph 11 MeV 
is seen. In calculations with expanded volume (8Vo), the system shows a sharp peak at 
T m 10.5 MeV. This peak is, however, not as sharp as the one for the heavier system. 
In calculations on limiting temperature in the model of liquid-gas phase equilibrium, the 
influence of Coulomb forces has often been emphasized |3|j2^| in the instability of the system. 
In the present calculation, we see a relatively small effect on the transition temperature. 
With the Coulomb force switched off, the transition temperature is shifted up by~ 1 MeV for 
both the confinement volumes 4Vq and 8Vo and the matter density becomes more uniform. 
This transition temperature is somewhat lower compared to the critical temperature for 
asymmetric nuclear matter || with isospin asymmetry equal to that of the nucleus. 

To summarize, we have calculated the caloric curve and the specific heat for two systems 
in a self-consistent Thomas-Fermi theory at two volumes, namely at 4 and 8 times the 
normal nuclear volume. The specific heat C v shows a peaked structure possibly signalling a 
liquid-gas phase transition at a temperature of ~ 10 MeV which is lower than the calculated 
critical temperature for infinite nuclear matter but larger compared to the calculated limiting 
temperature for finite real nuclei [p73 ] . In simplistic model calculations pl| , it has been shown 



that the inclusion of correlations brings in features reminiscent of a phase transition in a 
system when no phase transition is evident in the usual mean field calculation; it would 
therefore be interesting to see whether fluctuations with two-body correlations bring down 
the phase transition temperature obtained in our TF calculation. 
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FIGURE CAPTIONS 



Fig. 1 The proton density profile for the system 150 Sm calculated at four tempera- 
tures in the volume V = S.OV . The dashed, dotted, dash-dot and full lines correspond to 
temperatures T = 5, 9, 9.5 and 10 MeV respectively. 

Fig. 2 The proton density profile for the system 150 Sm calculated at temperature T = 10 
MeV in two different volumes. The full and dashed lines correspond to calculations at 

V = 4.0Vo and V = 8.0V , respectively. 

Fig. 3 The temperature plotted as a function of excitation energy per particle (caloric 
curve) for the system 150 Sm. The dashed curve corresponds to calculations with volume 

V = A.OVo while the full curve corresponds to V = 8.0V . 

Fig. 4 The specific heat per particle plotted as a function of temperature for the system 
150 Sm. The dashed curve corresponds to calculations with volume V = A.0V while the full 
curve corresponds to V = 8.OV0. 

Fig. 5 Same as Figure 3 for the system 85 Kr. 

Fig. 6 Same as Figure 4 for the system 85 Kr. 
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